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ABSTRACT 

This article is intended to describe the significance of lineaments (surface and subsurface lineament) for 
groundwater exploration in part of Vindhyan fringe belt of Mirzapur district , Uttar Pradesh , India. Remote sensing is 
vital visualization tool for identification of geological feature , such as lineaments , land use pattern , geomorphologic 
features etc., and also applicable for groundwater investigation. The lineaments (surface and subsurface lineament) play 
a crucial role in groundwater investigation in hard rock area. The surface and subsurface lineaments are extracted from 
the LANDS AT-TM satellite imagery and grided free air anomaly (satellite gravity data), respectively. The correlation 
between surface and subsurface lineaments are performed by using the rose diagram and understand the orientations of 
traverse lines of geophysical survey. The integrated lineament and lineament density maps are depicted the appropriate 
locations were selected for detailed geophysical survey and also further detailed study. 
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INTRODUCTION 

Lineaments are one of the essential topographic features used in exploring for various natural resources 
such as groundwater, hydrocarbons, minerals, and geothermal energy, as well as in mapping hazard susceptibility 
from earthquakes and landslides (Guild 1974; Dix and Jackson 1981; Sibson 1986a, b; Boucher 1995; Rowan and 
Bowers 1995; Rowland and Sibson 2004; Masoud et al. 2007). In this regard, the characteristics of lineaments 
such as spatial extent, density, intersection and orientation have proved significant because they indicate zones and 
trends of high permeability and / or low pressure that may serve as pathways for migration and thus targets for 
increased reserve. The lineaments may represent faults that control the basin development and distribution of 
reservoirs (Warner 1997). 

Generally, lineaments are underlain by zones of localized weathering and increased permeability and 
porosity. Meanwhile, some researchers studied relationships between groundwater productivity and the number of 
lineament within specifically designated areas or lineament density rather than the lineament itself (Hardcastle 
1995). Therefore, the mapping of the lineaments closely related to the groundwater occurrence and yield is 
essential to groundwater surveys, development and managements. Many groundwater exploration projects had 
higher success rates when sites for drilling or detailed geophysical survey were guided by lineament mapping 
based on the satellite images (Teme and Oni 1991; Guftafsson 1993). On the basis of technology, the lineament 
can be classified into two types: surface lineament and subsurface lineament. The surface lineament can be 
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extracted from satellite imagery and subsurface lineament can be mapped out from geophysical techniques. 

The Bouguer anomaly is calculated to apply atmospheric correction, Bullard-B correction and complete detailed 
terrain correction (distant, intermediate and inner zone correction) from global free air grided data (satellite gravity data) 
using FA2BOUG (FORTRAN CODE) with the input parameters (Fullea et al. 2008). The extraction of subsurface 
lineaments from Bouguer anomaly is carried out using filtering technique (gradient of gravity data in various directions) 
with the help of 3 X 3 gradient operators. The integrated lineament and integrated lineament density maps are prepared by 
overlaying the thematic maps of surface & subsurface lineament map. The integrated lineament and integrated lineament 
density map of study area show the appropriate locations which are selected for geophysical survey and also for further 
detailed study. 

GEOGRAPHY AND GEOLOGY OF STUDY AREA 

The area chosen for the study is the part of Vindhyan fringe belt which is located in Mirzapur district of Uttar 
Pradesh (Figure 1). The study area partly falls in Survey of India toposheet (No. 63K/16, 63L/13, 630/4, 63P/1) at 
1:50,000 scales and lies between latitudes 24°57’54” and 25°06’50”N and longitudes 82°50’23” and 83°04’22”E. It is 
situated in the extreme north-eastern part of the district of Mirzapur, Uttar Pradesh covering parts of three blocks of 
Mirzapur district such as Narayanpur, Jamalpur and Rajgarh. The district Varanasi lies to its north, district Chandauli to its 
south-east, and remaining part in the district Mirzapur. 

The area chosen for the investigation is a part of northern portion of Vindhyan fringe belt which is situated in the 
south-eastern direction of the river Ganges. Some part of the area falls on the northern slope of Vindhyan plateau which is 
entirely covered with a varying thickness of Pleistocene and Recent alluvial deposits of the Gangetic river system 
(Krishnan and Swaminathan 1959). The Vindhyan Supergroup is composed mostly of low dipping formations of 
sandstone, shale and carbonate, with a few conglomerate and volcaniclastic beds separated by a major regional and several 
local unconformities. The regional unconformity occurs at the base of the Kaimur Group and divides the sequence into two 
units: the Lower Vindhyans (Semri Group) and the Upper Vindhyans (Kaimur, Rewa and Bhander Groups). The outcrop 
pattern of the Supergroup resembles a simple saucer-shaped syncline. It is generally believed that the Vindhyan basin was 
a vast intra-cratonic basin formed in response to intraplate stresses (Bhattacharyya 1996; Bose et al. 2001). 

Groundwater in Kaimur Formation may occur in the weathered and fractured sandstone provided the zone is 
connected with recharging sources. The Kaimur Group of the Vindhyan Supergroup is of special significance because it 
consists dominantly of silica-clastic rocks lying unconformably over the carbonate-rich Semri Group-Lower Vindhyans. 
Therefore, the rocks from the Kaimur Group hold strong evidence regarding changing environment of deposition, climatic 
conditions & tectonics and weathering conditions, during Mesoproterozoic (Mishra and Sen 2010). 

The southern portion of study area is covered by thin surface soil with varying thickness ranging from 1-4 m and 
few exposure of Vindhyan sandstone whereas thickness of overburden increases while moving towards northern side. The 
bedrock (Vindhyan Supergroup of rocks) is expected below the surface soil cover (Krishnan 1982). 

MATERIALS AND METHODS 

Data Collection 

The LANDS AT-TM (Path- 142/ Row-043) data is used in the present study. The multispectral satellite image from 
the LANDSAT-TM sensor contains seven bands, i.e. Blue (0.45-0.515 pm), Green (0.525-0.605 pm), Red (0.63-0.69 pm), 
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Near Infrared (0.75-0.90 jam), Shortwave IR-1 (1.55-1.75 jam), Thermal IR (10.4-12.5 jam) and Shortwave IR-2 (2.09-2.35) 
(Palamuleni 2009). The survey of India toposheets 63K/16, 63L/13, 630/4 and 63P/1 with a scale of 1:50,000 are used as a 
source of ancillary information. 

The publicly available global grided free air gravity anomaly (satellite gravity data) and elevation data sets are 
used in the present study. Elevation data (topography) is taken from the 1 ' X 1 ' new release of Etopo2 2 global data base 
(V9.1) by Hastings and Dunbar (1998); Sandwell and Smith (1997) and Smith and Sandwell (1994). The grided free air 
gravity anomaly data is taken from the V X V global satellite altimetry data compilation VI 6.1 2 model by Sandwell and 
Smith (1997). 

Surface Lineaments Extraction 

Arkoprovo et al. (2012) discussed about the satellite data analysis, the main task in this stage is to carry out 
analysis and interpretation of satellite data, in order to produce thematic maps, such as lithology, structural and land use 
maps. Initially, all the images are rectified using the Survey of India (SOI) toposheet. This is followed by processing the 
digital images using the various processing techniques, viz., enhancement, filtering, classification and other GIS processes 
16. Subsequently, selective field checking is carried out. 

Subsurface Lineaments Extraction 

Fullea et al. (2008) explained to perform the terrain correction (Bullard C) in three steps: the distant zone 
computation, intermediate zone computation and inner zone computation. In a distant zone computation, the harmonic 
spherical expansion of the potential of each right rectangular prism representing an elevation grid point has been 
considered by Fullea et al. (2008). In an intermediate zone computation, Fullea et al. (2008) computed the gravitational 
attraction produced by each prism using the analytic formula. In case inner zone computation, Fullea et al. (2008) divided 
the scheme of the computation of inner zone contribution into two parts: (i) a flat- topped prism with a height equal to the 
elevation of the computation point; (ii) four quadrants of a conic prism sloping continuously from each square of the inner 
zone to the calculation point. Fullea et al. (2008) developed a computer program written in FORTRAN 90 (FA2BOUG) 
specifically designed to determine the Bouguer anomaly from publicly available global grided free air anomaly and 
elevation database sets. FA2BOUG computes the complete Bouguer correction (i.e. Bullard A, B and C corrections) in 
several spatial domains according to the distance between the topography and the calculation point. In the present study, 
FA2BOUG computer program has been used with input parameters to obtain a complete Bouguer anomaly map of the 
study area. Finally, mapping of subsurface lineaments are carried out using gradient of gravity data in north-south and east- 
west direction. 

RESULTS AND DISCUSSIONS 

Surface Lineament and Surface Lineament Density Map 

The satellite images due to its many capabilities such as the synoptic aerial coverage, multispectral captivity of 
data, temporal resolution, etc., produce better information than conventional aerial photographs (Lillesand and Kiefer 
1999), so the same has been selected for the task of extracting surface lineament. Radiometric correction is applied to 
images to remove noise present in the image. A 3*3 median filter is applied to all bands of LANDS AT-TM (Path- 142/ 
Row-043) satellite images and maintain the output values which are not affected by the actual values of outlier cells within 
the filter window. The median filter is particularly good for removing isolated random noise. With the help of image 
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interpretation keys and elements, lineaments are interpreted for the entire study area using standard False Color Composite 
(FCC). The digital image enhancement techniques can be contributed significantly in extracting surface lineament. Among 
the variety of image enhancement techniques, the filtering operations (Suzan and Topark 1998; Chang et al. 1998; Mah et 
al. 1995), Principal Component Analysis (PCA) (Qari 1991; Nama 2004) and spectral ratioing (Arlegui and Soriano 1998) 
are the most commonly used ones (Mohammed et al. 2010). The surface lineament map (Figure 2) has been created by 
using FCC and image enhancement techniques. Arkoprovo et al. 2012 described lineament density map. The lineament 
density map was prepared by dividing the study area into 1 km/ 1 km grids. The total length of lineament in each grid was 
measured and plotted in the respective grid centers. These values are joined by isolines to prepare a lineament density map 
as shown in Figure 3. 

Surface lineament density map was generated from the surface lineament map. The values of surface lineament 
density have been assigned depending on the density of the surface lineament. The high surface lineament density areas 
may have important of groundwater prospects, hence acts good groundwater prospect. 

Subsurface Lineament and Subsurface Lineament Density Map 

The extraction of subsurface lineaments from Bouguer gravity anomaly is carried out using gradient of gravity 

data along north-south and east-west direction with the help of 3 X 3 gradient operator. It is grided and contoured to 

produce gradient directional map of gravity data along north-south and east- west direction by Surfer Version 10 and 
Matlab Version 8.3. 

The gradients of gravity data along east- west direction are grided and contoured to produce a gradient directional 

map of gravity data along east- west direction. The subsurface lineaments are primarily oriented in one major set of 

direction i.e., approximately north-south direction. The maximum number of subsurface lineaments occurs in southern part 
of study area as shown in Figure 4. 

Further, the gradients of gravity data along north-south direction are grided and contoured to produce a gradient 
directional map of gravity data along north-south direction. The map shows that the subsurface lineaments occurs in 
southern part of the study area in more number as compared to those derived from gradient directional map of gravity data 
along east- west direction. The subsurface lineaments are primarily oriented in one major set of direction i.e., approximately 
east- west direction as shown in Figure 5. 

The subsurface lineament map is prepared by overlaying subsurface lineaments (as shown in Figure 4 & 5) from 
filtering technique (gradient directional filter of gravity data along north-south and east-west direction) as shown in Figure 
6. The subsurface lineament density map has been generated from the subsurface lineament map. The values of subsurface 
lineament density have been assigned depending on the density of the subsurface lineament which is shown in Figure 7. 
The high subsurface lineament density region may be identified as good prospects for groundwater resources. 

Integrated Lineament and Integrated Lineament Density Map 

The integrated lineament map is prepared by overlaying the thematic maps of surface and subsurface lineament 
map. The integrated map shows the appropriate locations which were selected for detailed geophysical survey and also 
further detailed study. It may or may not delineate the fracture zone because maximum numbers of lineaments are existed 
near the wells and tube wells as shown in Figure 8. Integrated lineament density map is generated from the integrated 
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lineament map. The values of integrated lineament density have been assigned depending on the density of the integrated 
lineament as shown in Figure 9. 

Orientation Analysis of Lineaments 

Linear feature orientation is considered as one of the most important characteristics. Rose diagram is essentially 
histogram for which the orientation axis is transformed into circle to give a true angular plot. The intervals of angle are 
plotted as pie-shaped segments of a circle in their true orientation and the length of the radius is proportional to the 
frequency of that orientation (Mohammed et al. 2010). 

Rose diagram is one of the most useful for representing the appropriate orientation of traverse line of geophysical 
survey and orientation of geological linear features. The surface lineament appears primarily in three sets of directions: two 
major sets oriented in northeast-southwest and northwest- southeast and the other one minor set oriented along 
approximately east- west direction as shown in Figure 10(A). The subsurface lineament appears primarily in two sets of 
directions: one major set oriented in east- west and other one minor set oriented along north-south as shown in Figure 
10(B). The integrated lineament appears primarily in four sets of directions: one major set oriented along east- west and 
three minor sets oriented in north-south, northeast- southwest and northwest-southeast as shown in Figure 10(C). It can be 
inferred from above study that traverse lines for geophysical survey have to be selected approximately in north-south 
direction depending upon the feasibility and availability of space in area. 

CONCLUSIONS 

The integrated lineament of study area depicts the appropriate locations which were selected for geophysical 
survey and also for further detailed study. It may or may not delineate the fracture zones because maximum numbers of 
lineaments or more density of lineaments are existed near the wells and tube wells. In integrated density map, the high 
density may have important groundwater prospects. The integrated lineament and integrated lineament density maps are 
very helpful for conducting geophysical surveys. These maps are very helpful in distinguishing the prospective sites to 
delineate hydrogeological features and the unfruitful attempts of conducting geophysical survey at several randomly 
selected locations could be avoided which significantly minimize the cost and time of geoelectrical resistivity survey. 
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Figure 1: Location Map of the Study Area 
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Figure 2: Surface Lineament Map from Satellite Imagery of the Study Area 
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Figure 3: Surface Lineament Density Map of the Study Area 
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Figure 4: Subsurface Lineaments from Gradient Directional Map of 
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Figure 6: Subsurface Lineament Map of the Study Area 
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Figure 7: Subsurface Lineament Density Map of the Study Area 
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Figure 8: Integrated Lineament Map of the Study Area 


www.tivrc.ors 


editor @tjprc. org 




Significance of Surface and Subsurface Lineaments for Groundwater Exploration 
in the Part of Vindhyan Fringe Belt of Mirzapur District , U. P., India 


23 



Figure 9: Integrated Lineament Density Map of the Study Area 
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